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Abstract: 

We present results for the light meson masses and decay constants as obtained from calcu- 
lations with the non-perturbatively improved ('Alpha') action and operators on a 24 3 x 64 
lattice at /3 — 6.2, in the quenched approximation. The analysis was performed in a way 
consistent with 0(a) improvement. We obtained: reasonable agreement with experiment 
for the hyperfme splitting; f K = 156 ± 17MeV, f v = 139 ± 22MeV, f K /U = 1-13(4); 
f K . = 219±7MeV , f p = 199±15MeV, = 235±4MeV; /J,(2GeV) = 178±10MeV, 
f p (2 GeV) = 165 ±11 MeV, where f v is the coupling of the tensor current to the vector 

mesons; the chiral condensate (qq) (2 GeV) = -(253 ± 25 MeV) 3 . Our results are com- 
pared to those obtained with the unimproved Wilson action. We also verified that the 
free-boson lattice dispersion relation describes our results very accurately for a large range 
of momenta. 



PACS numbers: 12.38.Gc,11.15.Ha,14.40.-n,13.30.E. 



1 Introduction 



From the very beginning of lattice QCD, one of the big challenges was to compute the 
hadron spectrum from first principles. In spite of the enormous technical progress that has 
been made in this field, yet there are ways to improve lattice studies systematically. The 
calculation of the light hadron spectrum is difficult mainly because of the large Compton 
wavelengths of the physical hadrons so that ever larger lattices are needed. On the other 
hand, to make a better contact with the continuum limit, simulations performed on several 
small lattice spacings are needed. These requirements are technically very demanding, and 
the search for systematic improvement is mandatory. 

Symanzik's proposal |l| for the improvement of the lattice action and quark bilinears 
with Wilson fermions, was realized perturbatively in B H]. Fairly recently, the Alpha col- 



laboration [4-10] (see also fllPp has proposed and to a large extent carried out, a thoroughly 
non-perturbative method which aims the elimination of all 0(a) discretization errors. In 
this way, one of the most important sources of systematic uncertainties in numerical studies 
on the lattice, is practically removed. 

The improvement program is implemented in several steps. The first source of O(a) 
uncertainties comes from the fermionic part of the Wilson lattice action. For on-shell 
quantities, these errors can be reduced to 0(a 2 ), by adding one higher-dimensional operator 
only. The resulting action reads: 

Ssw = Swuson + ic sw go ( y J2 Q( x ) cr ^F^(x)q(x) \ . (1.1) 

\ X,fJ,U / 

This is the Sheikholeslami-Wohlert or Clover action, where the last name is due to the 
shape of the lattice operator used for F^x). The non-perturbative determination of c sw , 
which is a function of the bare coupling only, allows the full non-perturbative improvement 
of the hadron spectrum. In Ref . || , c sw was determined non-perturbatively for different 
values of bare gauge coupling, and the final result of an overall fit for g$ < 1, is: 



1 - 0.656^ - Q.152^ 4 - 0.054^ 6 
Csw[%) ~ 1 - 0.922 gl • 



For (3 = 6.2 i.e. g$ = 6/(3 = 0.9677, this gives c sw = 1.614, which is the value used in the 
present studyQ. 

The second source of 0(a) errors, comes from discretization effects in the matrix ele- 
ments of composite local quark operators. These errors are relevant in the calculation of 
decay constants and/or form factors, i.e. quantities for which the knowledge of a hadronic 
"wave function" becomes crucial. As for the quark action, quark bilinears may also be 
improved through local counterterms, i.e. by adding specific operators which satisfy the 



1 In one-loop lattice perturbation theory at (3 — 6.2, to one-loop order c sw = 1.257, while with the so- 
called boosted coupling Jll], c sw = 1.479 (at tree level c sw = 1). With these c sw -values, calculations 
were already performed several times (a recent review with a complete list of references can be found in 



jj; see also (14|). 
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same symmetry properties as the original ones. Most of the coefficients of the countert- 
erms, which will be needed in this study, were calculated non-perturbatively. Their values 
will be given in course of the presentation. 

Some of the counterterms are present only out of the chiral limit and depend explicitly 
on the quark masses. They come with the so called ^-coefficients. The only one which is 
easy to obtain, from the forward matrix element of the vector current, is by, which has 
already been computed non-perturbatively |J. In order to determine ^-coefficients for the 
other operators, the improvement program was extended in Ref. and most recently 
in [0, but these proposals have not been applied yet to the "Alpha" action. Another 
attempt has been tried in Ref. [17 but the results are not stable [Q. For these reasons, 
we have taken the values of all the other 6-coefficients from perturbation theory. Their 
values, as well as those of the renormalization constants, will be quoted whenever used. 

The implementation of the improvement program in practical calculation was already 
done in [HJ and in [pOfl . As for the heavy quark sector, only preliminary numbers have 
been reported so far [21], and the final results will appear soon P2| . The main results of 
the present study, which concerns the light hadrons only, are given in the conclusion. 

This paper is organized as follows: in Sect. 2 we give a short outline of the lattice 
setup, compute the spectrum of light mesons and extract n crit in a way consistent with the 
improvement; in Sect. 3 we discuss the hyperfine splitting and the J-parameter; Sect. 4 is 
devoted to the study of decay constants; in Sect. 4, we also give our estimate for the chiral 
condensate; in Sect. 5, we make a comparison with previous (unimproved) lattice results; 
in Sec. 6 we test the energy-momentum relation on the lattice; we conclude in Sect. 7. 



2 Hadron Masses, Kcritic an d a 

In this section, we will briefly discuss the standard procedure for extracting hadron masses 
and fix the lattice parameters (^critic and a -1 ). We will insist on details only when the 
procedure is different as compared to previous (standard) analyses. 



2.1 Lattice Setup and Hadronic Masses 

Our results are based on a simulation performed on two Torre-APEIOO (25 Gflops) ma- 
chines located at the "Roma - I" University. Altogether, we have produced 100 gauge field 
configurations on a lattice of size 24 3 x 64 at (3 = 6.2, in the quenched approximation. 
After 5000 Metropolis sweeps, obtained by starting from a cold configuration, independent 
configurations were generated with a separation of 2000 sweeps. The values of the light 



Wilson hopping parameters used in our simulations, which are the same as in Ref. fL9 |, are 
the following ones: 

• 0.1352 (u); 0.1349 (d); 0.1344 (s); 0.1333 (1). 

The label assigned to the different quark masses (hopping parameters) are not to be 
confused with actual quark masses. The quark propagators were inverted using the minimal 
residual algorithm preconditioned a la Oyanagi p3[. 
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Figure 1: Effective masses for pseudoscalar (left - ) and vector (right figure) mesons. From 
up to down, four curves in each figure correspond to mesons containing 'IV, 'ss', 'dd', 'uu' 
quark flavours, respectively. 



To estimate the statistical errors, the raw results for various correlators (see below), 
as obtained on our 100 configurations, were jackknifed by decimating five of them at the 
timefl. 

Enabled by our correlators, we proceed the analysis by plotting the effective masses 
for pseudoscalar and vector mesons. These plots are shown in Fig. [I], where the effective 
masses are obtained as solutions of: 



coshM; 



tanhMr tanh 



T 

MA-r-t 



[2.1] 



The hadronic masses in lattice units will be denoted by Mj, whereas the physical ones by 
mj, i.e. Mj = amj. Time distances and space coordinates are always expressed in lattice 
units. By C JJ (t), we generically refer to the usual two — point correlation function summed 
over x, which (for large euclidean times) is dominated by the lightest hadronic state which 
couples to the chosen interpolating current J: 



CUi) = £(0|J(£,t)Jt(0)|0> 



t>0 



Zj 
2Mj 



-Mjt _|_ ^ e -Mj(T-t) 



^e- M ^ 2 cosh 

Mj 



Mj (? -t 



(2.2) 



T = 64 is the lattice temporal extension, and n the temporal inversion (t <-» T — t) sym- 
metry factor, which is +1 in the JJ-case for mesons^ Note that, this ('cosh') form has 

2 We also tested that by varying the number of configurations per cluster, the error estimates remain 
stable. 

3 Among the correlators considered in this study, only in the case C AP (t) (corresponding to the correlator 
of the fourth component of the axial current with the pseudoscalar density), one has: i] = (—1). 
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been used to obtain the relation Q2.1|) . For pseudoscalar mesons, the standard interpo- 
lating current that couples to the pion (J PC = CT + ) is Jps(x) = iq(x)j 5 q(x) . For the 
extraction of quantities related to the p-meson (J PC = 1 ), the local vector current, 
Jfi{. x ) — q(. x )lfi<l( x )i is the appropriate one. More specifically, we consider the space com- 
ponent Ji(x) and average over the indices, which is the procedure usually employed to 
reduce the statistical noise. By using effective mass plots, we may fix the initial time (ti n ) 
of the range on which we fit the data to extract the lightest masses^. The final time is best 
fixed by direct inspection of the signal to noise ratio in the hadronic propagator C 7J (t). 
With these two criteria, we establish the fit intervals: our light pseudoscalar mesons are 
well isolated for t £ [14, 29] , while the vector ones for t G [14, 24] . 

From the fit ( |2.2| ), we obtain the hadronic masses in lattice units Mj. The fit parameters 
are reported in Tab. [I], where we also present the results of our extrapolation to K crit which 
we discuss now. 



"flavor" 


M PS 


Zps 


My 




11 


0.4167(15) 


0.0111(4) 


0.4911(29) 


0.0037(2) 


ss 


0.3058(19) 


0.0077(4) 


0.4055(47) 


0.0022(2) 


dd 


0.2440(21) 


0.0063(4) 


0.3626(78) 


0.0016(2) 


uu 


0.2007(26) 


0.0057(4) 


0.335(12) 


0.0013(3) 


critical 




0.0035(4) 


0.275(22) 


0.0005(3) 



Table 1: Masses and Z 's for pseudoscalar and vector mesons in lattice units. These results 
are in good agreement with results of Ref. [l^ . 



2.2 Critical Parameter and Inverse Lattice Spacing 

We now discuss the uncertainties in the determination of K cr a which depend on the method 
used to fit the pseudoscalar meson masses. K crit represents the value at which the chiral 
symmetry on the lattice should be restored. In practical calculations, the basic relation is 



the Gell-Mann-Oakes-Renner one 24 



ml = -^(qq)m q + 0{m\) , (2.3) 

which states that the terms responsible for explicit chiral symmetry breaking are linear in 
the quark masses. On the lattice, for degenerate quark masses and by neglecting the terms 
of order 0{m 2 q ), this implies: 

M PS = arn q = | ( 1 - -L ) , (2.4) 



At that ti n , we assume that contributions of higher excitations which couple to a given correlation 
function are negligible. 
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where the standard definition (that can be derived from the vector Ward identity) of m q 
has been employed: 

am q = ~ (1 - J-V (2.5) 

However, this (standard) procedure to determine K cr u is valid if we have points sufficiently 
close to the chiral limit so that, up to chiral logarithms, higher-order quark mass terms 
can be neglected. These terms can arise from two sources. On the one hand, they are due 
to the lattice artifacts and can be eliminated by replacing am q —>■ am q (l + b m am q ). On 
the other, they can be a real physical effect, as indicated in ( |2.3| ). In order to investigate 
this point, we made a fit of the form: 

M PS = aA—- —\ + a 2 (— - —] (2.6) 

\ ™q ^crit J \ ^crit / 

from which we obtained our best estimate: 

• ^ = 0.135845(25) (2.7) 

with ol\ = 1.106(32), a 2 = 0.94(13). From the result of the fit, and as can be seen in Fig. |2], 
the sign of a 2 is opposite to what one would expect from the present determination of b m . 
For example, in lattice (boosted) perturbation theory b m = —0.593 (b m = —0.652) M []. 
Thus, unless perturbation theory gives the opposite sign (which we believe it is impossible), 
this implies that the positive curvature is a physical effect. The value of n crit Q2.7|) is the 
one that we will use throughout this and our forthcoming studies. We note, in passing, 
that the result of a linear fit with the three lighter mesons gives K cr u = 0.135801(19). An 
important observation is that ( |2.7|) agrees with the value K cr i t obtained from the axial 
Ward identity |25|: 

Pw < ~ 2 (P(t)Om) ( } ' ( ) 

In the above formula, the axial current was improved as A^(x) —> A^(x) + ca8^P(x) . A hat 
denotes that the quantity is properly renormalized (see Sect. 4). With the simple choice of 
0(t) = P(t), and by using the symmetric definition of the lattice derivative, at large time 
distances, we practically fit the following: 

Pwi — P (0) + c A a P (1) where 

(-4 ( f )P(0)) 
H 2(P(t)P(0)) V b! 

ap ( 1 ) = cosh (M PS ) - 1 (2.9) 



Results are presented in Tab. 2.2 



5 A non-perturbative estimate, b m = —0.62(3), was given in p7| . We tried to use the same technique 
but the values that we obtain are very unstable. 
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Figure 2: Quadratic vs. linear extrapolation of M ps to n crit . 





p (0) 


-c A apW/p (0) 


Pwi 


11 


0.0794(15) 


0.0410(8) 


0.0764(15) 


SS 


0.0450(15) 


0.0388(12) 


0.0432(15) 


dd 


0.0292(16) 


0.0378(19) 


0.0281(16) 


uu 


0.0199(16) 


0.0373(29) 


0.0192(17) 



Table 2: Lattice axial WI quark bare masses p's from mesons consisted of degenerate 'fla- 
vors'. The results are given in lattice units. Note that the 0(a) correction c A apW / 'p(°> 
never exceeds ~ 4%. 



The resulting value for n crit extracted from the linear dependence, p(l/n), is[J 

0.135840(48), (2.10) 



crit 



in very good agreement with Q2.7|) . Note also that ( p.lOj ) agrees with n crit = 0.135828(5), 
obtained by the same method in |20|| , as well as with the result of Ref. f[9 |: K cr n = 
0.13589(2). 

Evidently, these estimates differ sensibly (as expected), when compared to the (boosted) 



3 We checked that the result of the quadratic fit is indistingishable from the value we quote, k c 
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perturbative value, at /? = 6.2: K^f t T = 0.1374 [|, 0. 



We now discuss the calibration of the inverse lattice spacing. A conventional method 
to set the scale is obtained by extrapolating the vector meson mass to the chiral limit 
and compare the obtained result in lattice units, M p (see Tab.l.), to the experimental 
value m p = 768 MeV. Using a quadratic fit of My in the quark masses, we obtain: 
a -1 = 2.69(14) GeV . This value is in good agreement with our preferred value for 
a -1 = 2.75(17) GeV that we discuss in the next subsection. In the real world (unquenched 
QCD), to fix the physical scale, we can use any hadronic quantity e.g. f^,Tn p , ...[] In the 
quenched approximation however, different physical quantities can lead to different cali- 
brations of the lattice spacing. The calibration of a -1 from different quantities, (m p , tuk*, 
fn, fp etc.), differ by less than the quoted statistical errors. For this reason we are unable 
to study this systematic effect. 



2.3 'Lattice physical planes' - procedure 

The uncertainty due to the extrapolation to K cr u can be circumvent if we adopt the so called 



"method of lattice physical planes", which was proposed and used in Ref. [27]. One 
starts with a definition of two physical lattice planes, i.e. expresses My and some other 
physical quantity obtained by lattice calculations, generically denoted by $, as functions 
of Mp S . In the plane (M V ,M PS ), one looks for the point where a fit to the lattice data 



meets the curve My = C s ^yMj> 5 , with C s i fixed by its experimental value, m^/rriK- The 
point where the two curves cross determines Mk and Mk*- At the same value of M K in 
the plane ($,Mj, 5 ), one reads off the corresponding &k,k*, in the kaon sector. Similarly, 
by fixing Cm = m p /m n , we can obtain $ n , P , without direct extrapolation to K crit . In this 
study, all our data are fitted quadratically in Mp s . 

To fix the value of the inverse lattice spacing, one proceeds as follows. From the cross 



in the first plane (My, Mp s ), of fitted data with the curves My = C^J Mp S (with C = C s i 
or C = Cee), we get: 

M n = 0.0491(42), and M p = 0.279(24); 

M K = 0.180(12), and M K * = 0.321(21). (2.11) 
Using M p , we fix the scale as a -1 = m p /M p and obtain: a _1 (m p ) = 2.75(22) GeV . In 



27], it was argued that due to the fact that the mass of K* is in range of masses directly 
accessible on the lattice, m^* is the most suitable quantity for the scale fixing. We adopt 
this argument and get: 

cr l {m K *) = 2.75(17) GeV, (2.12) 

which is the value that we will use in the following discussion and all our forthcoming 
studies. 



7 In that case we would not use m p , since the p-meson is not a stable physical particle. 
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By using the same method, one can also determine the value of the light quark mass. 
Details of this analysis were presented in Ref. p6[ . Here, we only give the value of K, str , 
i.e. the one which corresponds to the strange quark mass: 

• K str = 0.13482(12) . (2.13) 

By using this result in Eq. ^6| , we obtain the hypothetical pseudoscalar: M 2 _ s = 0.0646(75). 
Then, in the plane (My, Mp s ), we read-off: 

Afy = 0.369(18) -> m^ = 1.013(19) GeV. (2.14) 

where we used a _1 (m^*). 



3 Hyperfine Splitting 

One of the main problems of lattice studies with Wilson fermions is to reproduce the exper- 
imentally observed, approximately constant hyperfine splitting. Theoretically, one expects 
my — m 2 PS to be constant in the heavy quark limit only. Experimentally, one finds: 



[m 



D* 



jrf^exp) = QM Q Ge y2 ; ( m 2^ _ ^(exp) = q_ 485 Q e y2 



There is no theoretical reason to explain why the hyperfine splitting has almost the same 



value also for light mesons |29 



(mi — m 2 



,(exp) = o.57i QeV 2 , (m 2 - m 2 ) {exp) = 0.552 GeV 2 . 



The net effect of the Clover term in the improved action is to give an extra anomalous 
magnetic moment to quarks, which increases and flattens the lattice hyperfine splitting. 
The most convenient method to study the splitting is to extract My — Mps directly from 
the ratio: 

%P-(t) ~ e^ M v- M rs)\ (3.1) 

The results for the mass difference and the hyperfine splitting are listed in Tab. |3], and 
displayed in Fig. ||. In this figure, we also show the results obtained with unimproved 
Wilson fermions (at the same j3 = 6.2 and on the same volume) [p7] . For comparison 
with experiment, we translated the experimental my — mp s , into lattice units. To this 
purpose, we take a _1 (cr) = 2.72(3) GeV, extracted from the computation of string tension 
|3tJ. We observe that the data obtained using the Clover action describe much better the 



experimental hyperfine splitting. 

It is encouraging that the Clover term reduces the slope in Mp s which characterizes 
the Wilson data. Moreover, the values are much closer to the experimental ones. Our 
preliminary study of the heavy hadron spectrum [2^] , however, shows that for heavy-light 



meson, all our results for the hyperfine splitting are well below the experimental values. 
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"flavor" 


My - M PS 


M v - M 2 PS 


11 


0.0735(22) 


0.0668(22) 


ss 


0.1006(42) 


0.0716(34) 


dd 


0.1219(83) 


0.0740(59) 


uu 


0.1398(140) 


0.0749(90) 



Table 3: Mass splitting in lattice units. 
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Figure 3: Hyperfine splitting on the lattice. The experimental line of results is given in 
lattice units by assuming a^ 1 (a) = 2.72(3) GeV. 



A consistent comparison between results for several physical quantities obtained with 
improved and the Wilson actions, including the hyperfine splitting, will be reported in 
Sect. 5. 

A frequently used quantity to test the effects of quenching (or other systematic errors) 
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is the dimensionless, so-called J-parameter |25|, |3l| : 



l PS 



J = M *iMT- (3 ' 2) 



In the kaon-sector, results obtained with the Wilson action are always below the experi- 
mental value, j( exp ) = 0.48. By fitting the My linearly in Mp s , the value for J that we 
obtain is: 

J = 0.373(7), (3.3) 



which is to be compared to J = 0.36(2) and J = 0.34(5) [^|, obtained on the same lattice 
and at the same j3 = 6.2, but with Wilson and tree-level improved (c sw = 1) Clover action, 
respectively. When we fit quadratically: My = AMp s + BMp s + C, we get: 



J = M K * (B + 2AM l K ) = 0.47(6) . (3.4) 

with A = -1.91(98)0, B = 1.57(31), and C = 0.275(22). Of course, the present statistical 
uncertainties, as well as the small number of hopping parameters that we use in this 
study, do not allow for a definite conclusion on this issue. In particular, the mesons with 
nondegenerate flavors are important to establish better the quadratic coefficient A, since 
the terms proportional to the square of the difference of the quark masses may give some 
contribution as well. Some more research to further investigate this point, is needed. 



4 Decay constants 

With the use of non-perturbatively renormalized improved operators, /„- becomes an equally 
good candidate as the hadron masses for fixing the inverse lattice spacing - a~ l . 

The standard procedure to extract the pseudoscalar and vector decay constants consists 
in calculating the following quantities: 

CpM = M3,t)P{0)) ^ p Jfo / T_ A 

C PP (t) Y*(P&t)P(0)) VZ^s V l 2 >) y } 

Cvv(t) = £(^(f,t)VK0)> ~ MyF v e~ Mv ^ cosh (m v £ - t)) (4.2) 

where the 'cosh'-form of fit ( |2.2| ) is assumed^. Note that we consistently use capital letters 
to emphasize that the quantity is given in lattice units; to distinguish it from the corre- 
sponding counterpart in physical units, which is denoted by lower-case letters. Fi refer to 

8 Note a small difference of A and the value reported in p6| , which is due to the shorter time interval 
for the fit of Cyy that we use in this study. This difference is however irrelevant for our final results. 
9 In the case of the pseudoscalar decay constant, we could also use Caa{£)- The reason why we use the 



one in Eq. (4.1), is that in this case the errors on decay constant are smaller. 
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P (o) 


a *PS 


Fps 






F v 


11 
SS 

dd 

uu 


0.0939(21) 
0.0832(27) 
0.0770(36) 
0.0733(47) 


0.1090(19) 
0.0891(21) 
0.0798(20) 
0.0754(21) 


0.0899(21) 
0.0799(27) 
0.0740(35) 
0.0706(46) 


0.1235(24) 
0.1161(32) 
0.1105(47) 
0.1061(69) 


0.0460(12) 
0.0349(14) 
0.0305(24) 
0.0303(46) 


0.1136(23) 
0.1086(32) 
0.1039(47) 
0.0996(70) 


critic 


0.0675(47) 


0.0636(23) 


0.0652(47) 


0.0996(88) 


0.0221(58) 


0.0949(87) 



Table 4: Decay constants - in lattice units. 



the decay constants of currents which include suitable overall renormalization constants. 
The pseudoscalar, F PS , and the vector, F v , decay constants are defined as usual: 

(0\A \PS(p = 0)) = iFpsMps and 

(0|V^(p = 0)> =iet ) F v My (4.3) 

The improvement of the axial and vector currents is achieved by adding the derivative of 
the pseudoscalar density, d^P, and the divergence of the tensor current, d u T^ u , respectively: 

(0\A \PS(p = 0)) - Z A (1 + b A am q ) [(0\A \PS) + c A (0\ad P\PS)} 

= tMpsFps = iM PS (F P 0) s + CAaPPs) (4.4) 
(0\V l \V(p = 0)) - Z v (1 + b v am q ) [(0\Vi\V) + c v (0\ad T i0 \V)} 

= iMyefPy = iM v ef ] {F§ ] + c v aF v l) ) (4.5) 

where for clarity, a in corrective terms on the r.h.s. are written explicitly, cy and ca are 
suitable coefficients, the values of which are chosen as to cancel 0(a) errors in physical 
matrix elements; the renormalization constant Za, h A , Zy and by will be discussed later 
on. By observing that: 

the corrections to the decay constants in actual calculations are then obtained using: 



aF$ = \- — sinh(M P5 ) and similarly, aF v L> = \- -sinh(My). (4.7) 
Mps My 

With all these relations, we extract the decay constants and list their values in Tab. [|. 
In that table, the improvement coefficient ca = —0.037 was used. It was obtained in |5j 
with an overall fit at g$ < 1: 

9 1 - 0.748 gl , , 

= - 0M ™* 1 - 0.977% - (48) 
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We also take cy = —0.214, as suggested by the preliminary study of p2| . Here, we did 
not account for the errors they quote, i.e. cy = —0.214(74). In the case of light quarks, 
the improvement term (aFy^) is very small anyway. Still, one comment is in order. Of 
all improvement coefficients, only cy differs by one order of magnitude from its (boosted) 
perturbative value cy = —0.026. For this reason, to be totally on the safe side, we have 
also calculated Fy with the perturbative cy given above. In this case the results for vector 
decay constants that we present in this paper are increased by about 5%. While this effect 
is less pronounced in the light hadron case (since the correction proportional to c v is rather 
small), it turns out to be very important in determination of the heavy- light vector meson 
decay constants [p2]| . 

We now discuss the values of the renormalization constants which have been used to 
obtain the physical currents and A^. 

The renormalization constants, Zy and Za, depend on the external quark masses, i.e. 
Z^y = Zji } y(m) . In the chiral limit, both constants were calculated non-perturbatively 
in ||, with the following results: 

1 - 0.7663 gl + 0.0488 ^ 

V ~ 1 - 0.6369 gl 

1 - 0.8496 gl + 0.0610 g* . . 

Za = 1 - 0.7332 so 2 ' ( 9) 

In our case, it gives Za = 0.8089, and Zy = 0.7927. The last step in relating a lattice 
decay constant to its continuum value is to account for the explicit quark mass corrections. 
The constant by has been computed non-perturbatively. The global fit, first given in ||, 
was updated in ||: 

1 - 0.7613 g 2 + 0.0012 g 4 - 0.1136 gt 
V 1 - 0.9145 gl 1 ' ' 

giving in our case: by = 1.404. The other constants are not known non-perturbatively 
and we have to rely on perturbation theory. To one loop accuracy, we have: bj = 1 + 
bfgl, with bfi = (7/0.11492(4), b^ = C/ 0.11414(4) @, = (7/0.10434(4)0, with 
Cf = (N 2 — 1)/2N. 

Boosted perturbative values are obtained by the replacement go — > gsQ- F° r comparison, 
in this way we get by = 1.242, rather close to the non-perturbative result. This makes 
us confident that by using for the constant 6^, the result of boosted perturbation theory, 
bA = 1.240, we committed negligible error for the light meson decay constants considered 
here. The corrective coefficients by t A enter with the bare quark mass am q defined in ( |2.5| ). 

In order to extract the desired physical quantities, we are supposed to make an ex- 
trapolation to the chiral limit as well as an interpolation to the strange mass sector. The 



results of a quadratic extrapolation to n crit are given in Tab. |[_]. This procedure is basi- 



10 g 2 B = gl/(P) where the average plaquette (P) = 0.6136 for /3 — 6.2 as inferred by our Monte Carlo. 

11 We checked that the extrapolation of decay constants remains the same regardless of whether we 
combine F^ ' 1 ^ before or after extrapolation. For heavy mesons, however, and must be combined 
before extrapolation |B2|. 
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Figure 4: Method of planes: In the plane a), dot-dashed line is obtained by fixing C s i while 
the dashed one by fixing Cm. Cross-points with the line which quadratically fits the data, fix 
Mk and M% at which one reads off the values of the decay constant in the plane b), where 
the quadratic fit has also been used. 



cally equivalent to the case of massless pion in the method of physical lattice planes. The 
method of lattice planes, as applied to our data, is illustrated in Fig. [|. 

• In the plane (F PS , Mp s ), we extract the pseudoscalar decay constants, using M n = 
0.0491(42) and Mk = 0.1801(115) as determined with the lattice planes method (see 
Eq. (|2.11|) ). The results are: 

F w = 0.0649(68) + c A 0.0653(25) = 0.0625(68) -> F n = 0.0506(55) ; 
F K = 0.0717(45) + c A 0.0732(23) = 0.0690(45) -> F K = 0.0568(36). 

We checked that the same result is obtained if we apply the method directly to Fp$. 
Converting these results in physical units, by using a _1 (m^.), we obtain our best 
estimate: 

f n = 139(22) MeV and f K = 156(17) MeV. (4.11) 

As a consistency check, we have also applied the method to the ratios Fps/My, and 
obtained: f w = 137(20) MeV and f K = 156(16) MeV. This kind of check is applied 
to other constants as well. However, it should be noted that we prefer to quote ( [4.1 1[ ) 
as our best estimate, because the lattice spacing is fixed uniquely (by uik*) for all 
the quantities considered. 

Note that, if we used the decay constants to calibrate the inverse lattice spacing we 
get a _1 (/ w ) = 2.59(29) GeV, and a~ l {f K ) = 2.82(18) GeV. 
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As for SU(3) breaking, from the direct ratio of the decay constants, we obtain: 



whereas by using 



fx/U- 1 = 0.123(55), (4.12) 



fK/mK * 1.06(4) (4.13) 



we get 



Our final estimate is then: 



U/rrif, 



f K /U~ 1 = 0.136(32). (4.14) 



- 1 = 0.13(4), (4.15) 



clearly below the experimental value, (fx/ f-w ~ l)* 6 ^ = 0.22. This result, neverthe- 
less, is bigger than the prediction of one-loop quenched chiral perturbation theory: 
f K /U- 1^0.07 



From the plane, (Fy,Mp s ), with the same criteria as for the pseudoscalar decay 

constants, we extract the vector decay constants: 

F p = 0.0724(96), and F K * = 0.0795(71), 

which in physical units, using as before a _1 (m^.), gives: 

f p = 199(15) MeV, and f K . = 218(7) MeV. (4.16) 

This is our best estimate. We have also applied the method to the ratios Fy/My, 

from which we obtain: 

f p = 199(11) MeV, and f K , = 218(7) MeV. 

We stress again that the reason to prefer the conservative result Q4.16D , is that the 
lattice spacing calibrated by rriK* is used for both F p and Fk*, which is clearly not 
the case for Fy/My. 

Similarly, from our data we predict 

'm 4 = 0.2310(52), i.e. U = 235(4) MeV, 



For comparison, the experimental values of decay constants are (see ]29| and |34|): 

f(exp) = 208 ( 2 ) Me y ; j{exp) = j{exp) = MgV ^ 

which are obtained from 1-prong r-decays (f p ± and fx*), and from electromagnetic 
decay widths ( f p o and fy). We see that in spite of the quenched approximation, our 
results are in very good agreement with experimental values. 

In the literature, one often encounters an alternative definition for the vector decay 
constant gy, which is related to ours as: gy = fy/My. For completeness, we give 
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the results in these units: 

g p = 0.260(14), g K , = 0.246(8) and = 0.231(6). 
These values agree with Ref. |35| . We close this subsection by some useful ratios: 



i = 1.43(31), tf- = 1.40(18), ^ = 1.10(5). 

Jn JK Jp 



4.1 Coupling with tensor current 

A phenomenologically relevant quantity is the tensor coupling to the vector meson ( J PC = 
1 ). The values of this coupling provide the normalization of the leading- twist wave 
function for transversely polarized vector mesons and play an important role in light- 



cone sum rule analyses for heavy to light (vector) meson decays. In Ref. |p6| , a value of 



fp(/ji) = 160(10) MeV (in MS scheme and at /i — 1 GeV), was obtained as best estimate, 



after considering different sum rule procedures intended to circumvent the problem of the 
pollution induced by the low lying positive parity (J PC = l +_ ) state, &i(1235), in the 
correlation functions relevant to the calculation. It is clearly desirable to have a lattice 
estimate of this quantity. We consider the same tensor-tensor correlator as in continuum, 
i.e. in ( |2.2|) we take J(x) = T pu (x) = q(x)cr pu q(x). With no-spatial momentum and by 
choosing T pu (x) — > T i0 (x) {i = 1,2,3), we project out only the negative parity states (in 
this way, no problem of mixing with positive parity states arises), and obtain: 

(0\f l0 \V(p= 0)) = le^F^My. (4.17) 

As in previous cases, the improvement program amounts to a redefinition 

T^{x) -> T M „(x) + c T {d^V v {x) - d v V^x)) 

where the new constant ct is known perturbatively only, ct = 0.00896(1)C/ g\ J8|. With 
our boosted coupling g B = 1.256, we have ct = 0.0188. Thus, the improved definition of 
( fTTTl ) reads: 

(0\f i0 \V(p = 0)) = ief ) M v Z T (n) (1 + b T am q ) (i^ (0) + ac T F y T(1) ) . (4.18) 

Analogously to the previous cases, out of chiral limit, we need bx, where we again rely on 
boosted perturbation theory, i.e. bx = 1.219 ||. 
In practice we have: 

if (0) = and aFfV = -^sinh(M y ). (4.19) 

The tensor current correlators were treated in the same fashion as vector ones, by aver- 
aging them over three Dirac indices. We also used the same criteria as before to fix the 
time-interval for fit, which in this case is slightly shorter (t G [13, 19]) than the one we used 
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a 


ss 


dd 


uu 


critic 


M T 


0.4914(25) 


0.4041(39) 


0.3616(58) 


0.3382(88) 


0.2833(97) 




0.0020(1) 


0.0011(1) 


0.0008(1) 


0.0006(1) 


0.0002(1) 



Table 5: Masses and Z 's of the vector mesons but with tensor current T^x) = q{x)a liv q{x) 
(in lattice units). 



for the vector correlators. For this case, the results of our fit for the mass and the bare 
coupling are listed in Tab. As expected, the results for M T are compatible with those 
for My, i.e. the masses obtained from the vector current correlators in Sect. 2 (see Tab. |I]). 
Using the method of lattice planes, for the tensor couplings we obtain: 



F = 0.0686(47) + c T 0.0267(62) = 0.0691(46), and 



F£. = 0.0721(30) + c T 0.0332(52) = 0.0728(30). (4.20) 

As for the renormalization constant, contrary to the axial and vector currents, 
is renormalization scheme and scale dependent. We remark, however, that the one-loop 
values of Zt in Landau gauge, in RI (MOM) and in MS schemes, are the same. We have 
computed the non-perturbative value of this constant in the chiral limit, in the RI (MOM) 
scheme, using the method of Ref. [37|. More details will be presented in [|18|| . Fig. |5], 



shows that for quark virtualities between 1 < (fia) 2 < 2, the one — loop evolution : 

/ a (a) \ 7 o T)/26 ° 

ZtM = ZtM (4 - 2i) 

describes very well the dependence Zt(h) on \x. The anomalous dimension of this operator 
is 7 = 2Cf, whereas b = 11 in the quenched approximation. We use ( |4 . 2 1| ) to get our 



non-perturbative value for Zt at /i = 2 GeV in the MOM (or MS) scheme: 

Z^ OM (/i = 2 GeV) = Z^Qjt = 2 GeV) = 0.87(2), (4.22) 

where the error is mainly statistical, plus the small uncertainty in inverse lattice spacing. 
This result is somewhat lower than the one obtained in (boosted) perturbation theory, 
Z T (2Ge~V) = 0.934(5) to one-loop accuracy g |9[ (with c sw = 1.614). 
Eventually, by using the relation Fy(fi) = Z|? s (/i)(l + bTam g )Fy and Eq. (|4.20|) , we get 
the following results: 

fj{2 GeV) = 165(11) MeV, and (2 GeV) = 178(10) MeV, (4.23) 
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1.1 



■ Non-perturbative Z T (u.a) 
— 1-loop evolution 




0.5 1.0 1.5 2.0 

(|ia) 2 

Figure 5: Zjn(fia) non-perturbatively computed on the same 100 configurations and with the 
same k's - extrapolated to the chiral limit which. Note that the one-loop evolution describes 
well the dependence Z T {\ia) for 2.7 GeV < /i < 3.9 GeV. 



where we took into account the statistical error and the error induced by the renormal- 
ization constant. This is our best estimate. As in previous cases, we made a consistency 
check and calculated these constants from Fy(fj,)/MT, and obtained: 

/ p T (2GeV) = 161(8), and /J*(2GeV) = 178(8) MeV. 

Other interesting quantities are the ratios of couplings of 'different' sources of vector 
mesons: 

(Ipj = 0.83(13), U^J = 0.82(8). (4.24) 



4.2 Chiral condensate 

Since we have all necessary ingredients, we can make an estimate of the value of the chiral 
condensate. A recent detailed discussion about different ways to extract this quantity from 
the lattice calculations, can be found in Ref. [[KJ . 

First, we rely on the GMOR relation (|2.3j ), where we use the quark mass defined by 
the vector ( |2.5| ) and/or the axial ( |2.8| ) Ward identities. The quark mass (|2.5|), as derived 
by using the vector Ward identity, was already used in the discussion on the determination 
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of K C rit- At this point, we can combine our result for F PS , with the coefficient from the 
fit (f2.6|), «i = 1.106(32). Then from ( |2.3| ), one easily obtains the relation 



{qq)^) = - l -Z s ^)a l Fl (4.25) 

which represents the first method (M-I) we use to extract the chiral condensate. F x de- 
notes the pseudoscalar decay constant, extrapolated to the chiral limit: F x = 0.0527(39). 
Equivalently, by taking the quark mass as derived from the axial Ward identity, we perform 
a fit analogous to the one in Eq. (|2.6|), i.e. 



M^ s = a 1 (2p)+« 2 (2p) 2 . (4.26) 

The values of the parameters are: <5i = 1.01(7), and a 2 = 0.81(35). This provides us the 
second method (M-II) to estimate the value of the quark condensate: 



m q(qq) is a renormalization group invariant quantity. Taken separately, the condensate 
is defined only in a specific renormalization scheme and at certain scale. These details 
are encoded in Zs,p{fi), which were recently calculated nonperturbatively [|]1| 26 1, in RI- 



scheme (MOM) and in Landau gauge. The results which we use here, are obtained after 
extrapolation to the chiral limit: Z s (2GeV) = 0.55(1), Z P (2GeV) = 0.43(1). With these 
values, we obtain the following results: 

(gg) (M) (2GeV) = -Z™ (jjl) (l.41(31) • 10~ 3 ) = -(244 ± 25MeV) 3 M-I 

4 R V) 



Z A 



1.30(21) ■ 10~ 3 = -(241 ± 24MeV) d M-II (4.28) 



where we used a 1 (m^») to express the result in physical units. 

For the conversion of our results from MOM to MS scheme, one relies on perturbation 



theory. At /i = 2 GeV, the matching was performed at the next-to-next-to-leading order 

f s NL ° = 1.243, and R N / ' 



in [p8|| . The value of the conversion factor is: Rg NL ° = 1.243, and Rg L ° = 1.144. Thus, in 



the MS-scheme, to NLO accuracy, we have: 



{qq) (2GeV) = -(255±26MeV) 3 VLO M-I 

= -(252±25MeV) 3 VLO M-II. (4.29) 

Another possibility (M-III) to estimate the value of the chiral condensate is provided 
by the low energy theorem: 

(0\qt l5 q\7t) = -A(0|gg|0>. (4.30) 
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For this purpose, we take the results for the matrix element |(0|P(0)|PS')| 2 = Zps, listed in 
Tab. (|I|). With its value extrapolated to the chiral limit, Z x , as well as with F x , Eq. (|4.30|) 
amounts to the following: 

(gg) (M) (2GeV) = -Zf ] {^^[z x = -(241 ± 20 MeV) 3 

-> (qq) MS (2 GeV) = -(252 ± 21 MeV) 3 LO M-III. 

Altogether, we quote 

(qqf S "(2 GeV) = -(253 ± 25 MeV) 3 LO , (4.31) 



as our final result^- Our estimates agree with results of Ref. 
Note also, that by using the evolution equation to NLO (with = 300 MeV), we obtain: 



(qqr (IGeV) = -(232 ± 24 MeV) 3 , which is very close to the value suggested by authors 



NLO ' 



of Ref. gl|, (qq) (1 GeV) ~ -(225 ± 25MeV)^ 

5 Comparison with unimproved results 



Action 


Wilson (c sw 


= 0) 


P 


6.0 6.2 


6.4 


# Confs 


320 250 


400 


Volume 


18 3 x 64 24 3 x 64 


24 3 x 64 


a- 1 (m p )[GeV] 
a-^mjrOtGeV] 


2.14(5) 2.84(11) 
2.26(5) 3.00(9) 


4.00(17) 
4.15(16) 



Table 6: Summary of lattice characteristics used for the comparison with the. a 1 {m p ) is 
obtained by the conventional - and a~ 1 (rriK*) by the lattice plane method. 



In this section, we make a 'consistent' comparison of our results with the unimproved 
data, which were generated by the APE collaboration in previous studies using the Wilson 
action. By 'consistent', we mean that the same methods to extract the physical quantities 
were used both in the improved and unimproved cases. In Tab. 0, we give some basic 



information on the simulation with the Wilson action. We refer the reader to Ref. |27| f° r 
more details. 



12 As indicated, all our estimates of (qq) S (2GeV) are given at NLO accuracy, so that it is easier to 
make a comparison with results of other approaches and other lattice groups. The results to NNLO can 
be obtained trivially, by multiplication by R s / R s 
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The physical volume for (3 = 6.0 and (3 = 6.2 is approximately constant ~ (1.7 fm) 3 , 
whereas for (3 = 6.4 it is somewhat smaller ~ (1.2 fm) 3 . Since all criteria used in Sec. 2. ,3. 
and 4. to extract masses and decay constants are applied to all lattice data, we can 
reliably investigate the effects of improvement. To extract physical observables, we employ 
the lattice plane method, since no explicit extrapolation to K, cr u is needed. The scale for 
each lattice is fixed by m K * and the results which are to be compared, are listed in Tab. |7|. 



Action 


Wilson 


Clover - NP 


a 


6.0 


6.2 


6.4 


6.2 


a(rriK*) [fm] 


0.087(2) 


0.066(2) 


0.048(2) 


0.072(3) 


U [GeV] 


0.160(7) 


0.138(7) 


0.150(10) 


0.139(22) 


fx [GeV] 


0.172(6) 


0.155(5) 


0.164(8) 


0.156(17) 


(Jk/U) - 1 


0.076(10) 


0.126(14) 


0.095(16) 


0.123(55) 


fp [GeV] 


0.310(14) 


0.286(16) 


0.241(12) 


0.199(15) 


fx* [GeV] 


0.312(10) 


0.291(10) 


0.255(8) 


0.218(7) 



Table 7: Comparison of several physical quantities calculated using Wilson fermions: non- 
improved and non-perturbatively improved. The lattice scale was fixed by ttik* ■ 



To have some better insight on the effect of the improvement, we plot in Figs. ^ four 
dimensionless ratios. We constructed these ratios using quantities directly extracted from 
correlation functions. From Figs. |6]a,b we see that the mass dependence of our results is 
very close to that obtained with the Wilson (unimproved) action at /3 — 6.4. However, this 
conclusion is only qualitative since the physical volume used for calculations at /3 — 6.4 
was small. 

As for the decay constants (Fig. |6|c,d), we use in both cases currents which are non- 
perturbatively renormalized. So the differences have to be attributed to genuine 0(a) 
effect. In the unimproved symptom of the presence of large 0(a) effects was 

particularly evident in the vector meson decay constant. Particularly problematic was 
the determination of Zy f4"2]] . Very recently in Ref. []4"3"| , Zy was calculated by using 
the program for non-perturbative renormalization |ft7|| , but the values for the vector decay 
constants remained well above the experimental ones. We compare our results in Figs. |6|c,d. 
We notice a significant change for the vector decay constant, less pronounced for the 
pseudoscalar one. This change improves the agreement of the quenched lattice results with 
the experimental values. We reiterate that the improved results would be about 5% bigger, 
had we used the value of cy from boosted perturbation theory. The reason to point this out 
is that the nonperturbatively determined Cy is almost by an order of magnitude bigger than 
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Figure 6: Comparison of improved (Clover-NP) with unimproved (Wilson) results, a) and 
particularly b) show the effect of improvement for the spectrum, c) and d) show this effect 
for decay constants. 



the corresponding perturbative value. While this observation does not substantially alter 
our results in the light meson sector of light meson, it seriously affects the determination 
of the heavy-light vector meson decay constants |22]. 

To better monitor the effect of the improvement, we also make the comparison of the 
vector decay constants with the results obtained by using the tree-level improved Clover 
action (c 3w = 1), at the same j3 = 6.2, and the same volume 24 3 x 64 (see [f27| for details). 
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In Fig. 7, we observe that the effect of the tree-level improvement is a slight decrease of 
fy. Further decrease towards the experimental values is the effect of the full elimination 
of 0(a) errors. 
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Figure 7: The effect of improvement on the vector meson decay constants. All results 
are obtained at f3 = 6.2 and the volume 24 3 x 64, with unimproved (c sw = 0), tree-level 
improved (c gw = 1), and non-perturbatively (c sw = 1.614) improved Wilson action which 
are denoted as Wilson, Clover-tree and Clover-NP, respectively. 



6 Energy momentum relation 

For studies of semileptonic and radiative decays on the lattice, it is important to calculate 
the form factors as functions of the momentum transfer. This is achieved by giving differ- 
ent momenta to the interacting hadrons. The injection of momentum introduces further 
discretization errors, in particular, affecting the continuum dispersion relation: 

E(Pf = P 2 + M 2 . (6.1) 

In order to investigate lattice artifacts related to this problem, we studied the meson 
propagators at several values of P. For this purpose, we used the Eq. (|6.1|) , as well as the 
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lattice free boson dispersion relation: 




sinh* I ^p- I = smh 2 ( — ) , Vshr I - I (6.2) 



which may be derived from the discretized free boson action with nearest neighbors inter- 
action. This is not a unique choice since it depends on the way we define derivatives on 
the lattice. In previous studies p7], H3, the authors have shown that the choice (|6.2|) is 



indeed favored by data. 

We used lattice cubic symmetry, parity and charge conjugation to relate symmetric 
configurations and thus increase the statistical quality of our correlators. By writing 
P = 2tt/ La (n x ,n y ,n z ), the components of spatial momenta for pseudoscalar and vector 
correlation functions that we consider here, are: 

(0,0,0); (1,0,0); (1,1,0); (1,1,1); 

(2,0,0); (2,1,0); (2,1,1); (2,2,0). (6.3) 

We limited ourselves to \P\ 2 < 8, after a test-run where we observed that for higher mo- 
menta, the hadronic propagators were immediately drowning into very large noise. We 
used the same criteria as in Sec. 2 to establish the time intervals for the fit at each momen- 
tum considered. These time intervals are listed in Tab. [S[ Note that for large momenta, 
the fit intervals become very short. 



— * 

\p\ 


1 


V2 


v/3 


2 




v/6 


VE 


PS ■ / — / 

1 kJ • u min v max 


10-23 


10-19 


10-18 


10-15 


10-14 


10-13 


10-12 


VEC . t m i n t max 


11-24 


11-21 


11-18 


11-18 


11-17 


11-15 


11-14 



Table 8: Values of the times for fit for pseudoscalar and vector correlation functions for 
all momenta considered in this work. 



To check the relations (|6.1|) and ( |6.2|) , we fit our data either with 



C JJ (t,P) = -^ s -e- E t p *, (6.4) 



or with 



Cjj(t, P) = ^h^ e ~ E{P)t - ( 6 - 5 ) 

2 sinh [E{P)j 

The same relations hold for transversely polarized vector mesons which we consider here. 
In general, however, the term due to the polarization should be accounted for. 
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Figure 8: Energy-momentum relation for light-pseudo scalar (a) and vector (b) mesons: In 
both figures E 2 vs. \P\ 2 as obtained from the fit to (|6.5|) , are denoted by bullets; [+] mark 
values obtained by QS.lp , while [□] denote values obtained by using ( |b.2p , where the masses 
are listed in Tab. [l|. 
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I p|2 _ n 

\r\ — U 


1 p|2 i 


1 PI 2 _ O 


1 P|2 _ Q 


vector 11 


0.4911(29) 


0.5576(24) 


0.6143(22) 


0.6680(21) 


pseudoscalar // 


0.4167(15) 


0.4957(18) 


0.5520(20) 


0.6190(23) 


vector ss 


0.4055(47) 


0.4839(36) 


0.5478(33) 


0.6099(33) 


pseudoscalar ss 


0.3058(19) 


0.4041(22) 


0.4829(24) 


0.5493(34) 


vector dd 


0.3626(78) 


0.4493(58) 


0.5168(49) 


0.5838(49) 


pseudoscalar dd 


0.2440(21) 


0.3558(27) 


0.4409(32) 


0.5077(49)) 


vector uu 


0.335(12) 


0.4281(74) 


0.4976(69) 


0.5678(67) 


pseudoscalar uu 


0.2007(26) 


0.3242(40) 


0.4156(41) 


0.4817(70) 




1 p|2 A 


1 P|2 C 


1 P|2 R 


1 p|2 o 


vector 11 


0.7084(22) 


0.7518(22) 


0.7933(33) 


0.8623(35) 


pseudoscalar // 


0.6678(28) 


0.7138(34) 


0.7556(48) 


0.828(7) 


vector ss 


0.6515(32) 


0.6987(33) 


0.7454(50) 


0.8165(45) 


pseudoscalar ss 


0.6029(46) 


0.6527(52) 


0.6989(83) 


0.779(15) 


vector dd 


0.6253(32) 


0.6729(44) 


0.7225(71) 


0.7931(56) 


pseudoscalar dd 


0.5720(79) 


0.6198(82) 


0.6712(139) 


0.767(31) 


vector uu 


0.6102(57) 


0.6572(55) 


0.7116(89) 


0.7783(67) 


pseudoscalar uu 


0.5633(152) 


0.6029(108) 


0.6595(220) 


0.779(65) 



Table 9: Values of energies E(P), as extracted from the fit to our data for different spatial 
momentum injections and for fixed Z = ZiP = 0)). 

On the basis of the Wigner-Eckart theorem, we fix Zj at \P\ — 0, to fit E(P). We 
checked and realized that the relation ( |6.5| ) provides a much better fit to our data, starting 
from the three-momenta \P\ 2 > 3. By using (|6.5|) , we find that the agreement with ( |6.2|) 
is excellent, up to rather large values of injected momenta. The results of these fits are 
given in Tab. ||, and illustrated in Fig. |[ In that figure, we depict our data by 'bullets' 
describing the results obtained from the fit to ( |6.5| ). For better comparison with dispersion 
relations, we also show E(P) as obtained from (|6.1| ) and ( |6.2|) and represent them by 'plus' 
and 'square' symbols respectively. The mass terms in these dispersion relation are those 
extracted from the correlators with \P\ = 0, i.e. those which we already listed in Tab. |l|. 
A comparison with previous results is somewhat difficult since the simulations with the 
Wilson action were done for \P\ 2 < 4 only. This makes it difficult to distinguish which 
dispersion relation is better to use. It is important to note that our data are also compatible 
with both dispersion relations, ( |6.1|) and ( |6.2|) , when the small momenta are considered. 
Only for large momenta it becomes clear that the Eq. ( |6.2|) describes our data much better. 
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7 Conclusion 



In conclusion, we summarize the main results of this analysis. We performed a lattice study 
of the light hadron spectrum and meson decay constants with non-perturbatively improved 
action and operators. We showed that, for directly accessible (not so light) quark masses, 
the physical contribution of quadratic quark mass terms to the mass of the pseudoscalar 
mesons, exceeds the effect of lattice artifacts. Thus we conclude that, in the quenched 
case, there is a positive contribution to m 2 PS in m 2 q . For all the quantities considered here, 
we determined the physical values using the method of physical lattice planes. Our main 
results are the following: 

- In spite of the use of the quenched approximation, we find that the values of inverse 
lattice spacing, as obtained from different physical quantities, are compatible with 



Ref. 



a- 1 (a) 



2.72(3) GeV 



this work 



a 1 {m K *) a 1 (m p ) a 1 (f K ) a 1 (f n ) 

2.75(17) GeV 2.75(22) GeV 2.82(18) GeV 2.59(29) GeV 



each other (within the errors), and with the one obtained from the string tension^. 
Since the least extrapolation is needed for Mk*, we decided to fix the scale by this 
quantity. 

- We verified that the hyperfme splitting in the region of light mesons is well reproduced 
by our data. We have also extracted J = 0.37, but noticed that by using the quadratic 
fit My(Mp S ), one gets J = 0.47(6) which is (despite its large error) closer to the 
experimental value. This point deserves further investigation. 

- For the pseudoscalar decay constants, we have: 

f K = 156 ± 17MeV, f w = 139 ± 22MeV, whereas the SU(3) breaking is (in spite of 
improvement) smaller than the experiment, i.e. (fx — fw)/ fw — 0.13(4). 

- For the vector decay constants, we obtain: 

f K . = 219 ± 7 MeV, f p = 199 ± 15 MeV, / = 235 ± 4 MeV. 

- In determination of the coupling of the vector meson to the tensor current, we also 
calculated Zt(/j) non-perturbatively, which is a new result. In Landau gauge, and at 
ji = 2 GeV, we extracted: 

Z lp=Zf OM = 0.87(2). 
Results for tensor couplings are: 

f K « (2 GeV) = 178 ± 10 MeV, f p (2 GeV) = 165 ± 11 MeV. 

- By using the Gell-Mann-Oakes-Renner formula, we estimated the value of the chiral 
condensate: 



(gg) MS (2GeV) = -(253 ± 25 MeV) 

13 It is worth mentioning that the inverse lattice spacing in Ref. 1 30 , was obtained by using y/a 
440 MeV. 
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- By using the same criteria of analysis, we consistently compared ours to the results 
obtained with unimproved Wilson fermions. We conclude the following: 

o Our meson masses and hyperfme splitting, as directly extracted from correlation 
functions, are qualitatively comparable to those obtained without improvement 
at $ = 6.4, 

o The contribution of 0(a) term in improved operator is of order of 5%. A striking 
effect of the full 0(a) non-perturbative improvement is evident in the case of 
the decay constant for vector mesons, and less so for the pseudoscalars. 

- We explored a wide range of momentum injections to study the energy-momentum 
relation for pseudoscalar and vector mesons. The lattice artifacts become important 
for higher momenta, but we show that the lattice dispersion relation for a free boson 
fits our data to excellent accuracy. This conclusion will be very useful for the study 
of semileptonic decays of heavy to light mesons. 
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